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1.  INTRODUCTION 

Most  inertial  navigation  systems  in  missile  and  projectile  applications  use  an 
integrated  approach  that  combines,  for  example,  location  data  derived  from  GPS  and 
orientation  and  acceleration  data  derived  from  onboard  sensors.  To  date,  the  use  of 
magnetometers  to  assist  with  the  determination  of  a  projectile’s  motion  or  heading  has 
been  severely  limited.  Two  obvious  reasons  are:  low  sensitivity  of  standard  compasses 
and  proximity  to  ferrous  objects  that  modify  the  strength  and  direction  of  the  earth's 
magnetic  field.  However,  this  is  not  the  case  for  free  flight  and  orbital  applications.  For 
example,  ARL(Aberdeen)  has  recently  utilized  a  non-MEMS  magnetometer  to  measure 
the  spin  of  artillery  projectiles  in  free  flight  [1],  Previously  such  data  had  been  collected 
using  an  optically-based  yawsonde  that  measures  angles  with  respect  to  the  sun  [2].  The 
magnetometer  sensor  utilized  a  bridge  circuit  comprised  of  giant  magnetoresistive 
(GMR)  elements.  When  the  bridge  (two  GMRs  shielded  and  two  unshielded)  sweeps 
past  a  field  line,  a  sinusoidal  output  is  generated.  The  duty  cycle  is  related  to  the  spin  rate 
of  the  object,  but  sensor  output  and  conditioning  electronics  are  limited  to  spin  rates 
above  50  Hz.  It  was  noted  during  these  flight  tests  that  amplitude  modulations  of  the 
GMR  signal  were  correlated  with  the  yawsonde  solar  angle  data.  Further  analysis  has 
shown  that  with  greater  sensitivity  and  the  efficient  placement  of  multiple  sensors  at 
different  orientations  on  the  same  device,  the  feasibility  of  high  accuracy  roll  attitude  and 
angular  orientation  measurement  is  extremely  high  [3].  A  multiple  sensor  device  is 
needed  to  avoid  nulls  that  will  occur  when  the  sensitive  axis  aligns  with  magnetic  field 
lines.  This  multi-sensor  requirement  leads  naturally  to  the  necessity  of  miniaturization 
and  MEMS-based  implementation.  Furthermore,  a  sensor  that  directly  measures  angles 
will  be  significantly  more  accurate  compared  to  angular  rate  sensors  that  are  presently 
used  for  many  inertial  measurement  applications.  The  accumulated  errors  (over  long 
flight  paths)  of  the  current  angular  rate  sensors  yield  position  and  attitude  estimates  that 
are  often  very  crude  and  that  rapidly  degrade  with  time. 


Fig  1.  Principle  of  operation  of  xylophone  bar  magnetometer 


Before  the  start  of  the  project,  JHU/APL  had  developed  a  novel  magnetometer, 
based  on  a  classical  resonating  xylophone  bar,  that  was  intrinsically  linear  and  had  a  wide 
dynamic  range  such  that  it  could  measure  magnetic  field  strengths  from  less  than 
nanoteslas  to  Teslas  [4].  Furthermore,  it  was  ideally  suited  to  miniaturization  using  a 
variety  of  MEMS-based  technologies.  This  should  result  in  magnetometer  arrays  with 
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miniature  sensor  elements,  highly  integrated  electronics,  inexpensive  cost,  and  well 
suited  to  meet  the  intended  application.  The  resonating  xylophone  bar  magnetometer 
uses  the  Lorentz  force  as  a  measure  of  the  magnetic  field  vector.  It  offers  the  potential  of 
being  light-weight  and  operating  in  a  low  power  mode.  True  to  its  name,  the  device  is 
based  on  a  classical  resonating  xylophone  bar.  The  principle  of  operation  is  shown  in 
figure  1.  The  active  element  is  a  conductive  bar  supported  at  the  nodes  of  its 
fundamental  mode,  of  frequency  fo,  of  mechanical  vibration.  An  alternating  sinusoidal 
current  [Icos(2jrft)]  is  fed  along  the  bar  through  the  supports.  In  the  absence  of  a 
magnetic  field  the  bar  will  remain  stationary. 

In  the  presence  of  a  magnetic  field,  B,  the  Lorentz  force  F,  normal  to  the  bar  and 
the  drive  current,  is  given  by: 

F  =  /xf?cos(27t/0  (1) 

This  force  causes  the  xylophone  bar  to  vibrate,  with  a  displacement,  d,  at  the  mid¬ 
point  of  the  bar  being  given  by: 


(2) 


The  static  deflection  of  a  beam  simply  supported  at  both  ends  due  to  the  Lorentz 
force  uniformly  distributed  along  its  length  is  given  by  [5]: 


d* 


5  FI  2 
384 EIb 


(3) 


where  h  =  ab3/ 12. 

The  phase  angle  between  the  driving  current  and  the  deflection  is  given  by: 
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arctan 
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In  the  case  of  a  free-free  bar,  where  the  effect  of  the  supports  can  be  neglected, 
the  resonance  frequency  of  the  fundamental  mode  of  the  bar  is  given  by: 


_  22.4  I  EIb 
Jo  ~ 


2n  \p wtl4 


(5) 


It  can  be  seen  from  these  equations  that  at  resonance  (/= f0),  the  deflection  of  the 
xylophone  bar  is  linearly  proportional  to  the  amplitude  of  the  drive  current,  7,  the 
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mechanical  Q  of  the  structure,  and  the  magnetic  flux  density,  B.  Furthermore,  since  the 
xylophone  bar  includes  no  magnetic  material,  the  sensor  doe  not  saturate  in  the  presence  of 
large  magnetic  fluxes.  This  results  in  a  magnetometer  with  an  extremely  wide  linear 
dynamic  range.  Macro-devices,  chemically  milled  form  Cu-Be  foils,  with  dimensions  of  5 
mm  x  0.5  mm  x  0.025  mm,  Q  factors  in  the  range  1000  -  5000,  and  drive  currents  of  the 
order  of  one  amp,  have  demonstrated  dynamic  ranges  in  excess  of  80  dB  with  noise  floors 
below  0.1  nT/VHz  [6]. 

Two  approaches  to  implementing  the  MEMS-based  xylophone  magnetometer 
were  pursued  on  this  project.  The  first  used  polysilicon  devices  fabricated  on  the 
DARPA  supported  Standard  Multi-User  MEMS  Processes  (MUMPs)  foundry  service. 
The  sensitivity  of  these  devices  was  predicted  to  be  limited  by  the  current  carrying 
capacity  of  the  polysilicon.  Therefore,  the  second  approach  was  based  on  replacing  the 
polysilicon  with  a  platinum/lead  zirconate  titanate/platinum  (Pt/PZT/Pt)  sandwich 
structure  for  the  xylophone  bar.  The  PZT  material  deposition  and  processing  was 
undertaken  at  ARL,  Adelphi.  Shock  testing  was  undertaken  at  ARL,  Aberdeen. 
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2.  POLYSILICON  MAGNETOMETER  DESIGN  AND  CHARACTERIZATION 

Throughout  the  course  of  the  contract  several  designs  of  the  polysilicon  vibrating 
bar  magnetometer  were  fabricated  using  the  3-layer  foundry  service  at  JDS  Uniphase  (ex 
MCNC  and  Cronos).  Each  design  was  divided  into  4  quadrants,  with  each  quadrant 
containing  different  series  of  layouts.  Where  possible  (i.e.  vibrating  bars  containing 
polyl)  the  bars  had  small  dimples  placed  in  the  center  and  near  the  ends  to  alleviate 
stiction  problems.  The  support  electrodes  are  anchored  either  to  the  polyO  or  silicon 
nitride  using  a  “5  of  diamonds”  pattern  to  minimize  the  loss  of  vibrational  energy.  [7] 


A  schematic  diagram  of  a  typical  polyl  device  is  shown  in  figure  2.  The  lower 
electrode  consists  of  a  0.5  pm  thick  polyO  layer  patterned  on  the  silicon  nitride-coated 
silicon  substrate.  The  xylophone  bar,  support  electrodes,  and  mounting  pads  are 
fabricated  from  2  pm  thick  polyl  suspended  (after  release  of  the  sacrificial  silicon 
dioxide)  2  pm  above  the  nitride  layer  or  lower  electrode.  The  mounting  pads  are  attached 
to  the  lower  electrode  by  the  indicated  patterns  to  provide  a  rigid  anchor  which  minimizes 
vibrational  coupling  with  the  xylophone  bar.  Contact  between  the  xylophone  bar  and  the 
nitride  layer  is  minimized  by  having  0.75  pm  deep  “dimples”  in  the  center  and  toward  the 
ends  of  the  xylophone  bar.  An  SEM  micrograph  of  a  released  and  mounted  500  x  50  x  2 
pm  xylophone  bar  is  presented  in  figure  3.  The  micrograph  clearly  displays  the  anchor 
pattern,  the  etch  holes,  and  the  topographical  evidence  of  the  dimples.  It  also 
demonstrates  the  clean  release  of  the  xylophone  structure. 

Evaluation  was  carried  out  in  a  purpose-built  evacuated  chamber,  using  a  bench- 
top  beam  deflection  microscope.  The  beam  of  a  laser  diode  is  directed  through  a  lens 
onto  a  node  position  on  the  xylophone  bar  in  a  way  that  the  reflected  beam  is  focused  on 
the  position  sensitive  detector  (PSD).  In  this  arrangement,  the  PSD  is  sensitive  primarily 
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Figure  3.  An  SEM  micrograph  of  a  released  and  mounted  500  x  50x2  /Jm 

xylophone  bar 


to  rotation  of  the  beam  at  the  node  and  not  to  vertical  displacement.  With  the  optical  lever 
of  26  cm,  the  achieved  sensitivity  was  of  the  order  of  AR/R=0.06/mrad  or 
AR=0.3V/mrad,  where  a  typical  R=5V  is  the  amplifier  sum  output  for  the  reflected  beam 
and  AR  is  the  difference  output.  Magnitude  and  phase  data  were  obtained  in  a  mT  static 
magnetic  field  by  scanning  and  detecting  at  the  frequency  of  the  sinusoidal  drive  current 
using  an  EG&G  Model  7265  Lock-in  Amplifier. 

Results  obtained  by  scanning  the  frequency  of  the  current  through  a  500  x  2  pm 
xylophone  bar  with  4  pm  and  10  pm  wide  supports  and  at  a  pressure  of  35  mT  are  shown 
in  figure  4  a  and  b,  respectively.  The  resonance  frequencies  and  Q-factors  of  the  two 
devices  are  78.1  kHz  and  6700,  and  95.6  kHz  and  7500.  For  reference,  the  predicted 
frequency  for  a  polysilicon  free-free  resonating  bar  (assuming  a  Young’s  modulus  of  160 
GPa)  is  69.2  kHz.  The  higher  measured  values  were  attributed  to  the  torsional  resistance 
offered  by  the  support  beams  [8]. 


Figure  4.  Frequency  scans  for  a  500x2  pm  xylophone  bar  with  4  pm  (a)  and  10  pm 
wide  (b)  supports  and  at  a  pressure  of  35  mT. 


A  mathematical  model  of  the  xylophone  magnetometer  incorporating  an 
approximation  for  the  mass  and  stiffness  contributions  of  the  support  arm  was  developed 
to  test  the  validity  of  the  above  notion.  Essential  elements  of  the  model  are  shown  in 
figure  5.  The  two  support  beams  are  modeled  as  vertical  and  torsional  springs  (with 
spring  stiffness  Ky  and  Kj,  respectively)  and  a  concentrated  mass,  M. 


The  dashed  line  represents  the  fundamental  flexural  mode  of  vibration  of  the 
beam.  There  are  two  lower  frequency  modes  associated  with  the  flexural  vibration  of  the 
support  arm,  but  these  are  not  of  interest  in  this  study  as  there  is  negligible  deformation 
of  the  resonator  bar. 

In  the  current  study,  the  resonator  is  modeled  as  a  beam  based  upon  the  Bemoulli- 
Euler  theory.  Use  of  this  theory,  rather  than  the  more  complex  Timoshenko  theory,  is 
justified  because  the  wavelength  of  vibration  is  much  larger  than  the  thickness  of  the 
beam.  When  the  vibration  wavelength  is  of  the  same  order  as  the  beam  thickness,  the 
effects  of  shear  deformation  and  rotary  inertia  included  in  the  Timoshenko  theory  become 
important.  The  partial  differential  equation  governing  the  motion  of  a  Bemoulli-Euler 
beam  is  given  by 


d4w(x,0 

dx4 


+  pA 


d^wjxA) 
dt 2 


(6) 


where  w(x,t)  is  the  transverse  displacement  at  any  position  x  along  the  beam  and  at  any 
time  t,  I  is  the  moment  of  inertia,  and  A  is  the  cross-sectional  area  of  the  beam.  The 
modulus  of  elasticity,  E,  and  mass  density,  p,  for  polysilicon  are  160  GPa  and 
2320  kg/m3,  respectively.  The  beam  resonator  has  twenty  5  pm  x  5  pm  holes  spaced 
along  its  centerline  to  assist  in  the  release  of  the  polysilicon  structure.  The  mass  density 
of  the  polysilicon  was  reduced  by  2.6%  to  account  for  this  mass  loss,  but  the  effect  on  the 
stiffness  is  considered  negligible.  Solution  of  Eq.  (6)  requires  the  resonator  be  divided 
into  three  segments:  the  two  cantilevers  and  the  central  span  between  the  two  support 
arms.  This  leads  to  twelve  simultaneous  equations  governing  force  and  displacement 
compatibility  between  the  three  beam  segments  and  for  the  two  free  ends.  The 
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frequencies  that  cause  the  determinant  of  the  resulting  12  x  12  matrix  to  vanish  are  the 
natural  frequencies  of  the  system.  However,  because  the  two  low  frequency  support  arm 
modes  are  not  of  interest,  the  12  x  12  matrix  can  be  reduced  to  a  5  x  5  matrix  [see  Eq. 
(2)].  This  significant  reduction  in  matrix  size  is  accomplished  by  taking  advantage  of  the 
fact  that  the  support  points  are  at  the  nodes  of  the  fundamental  flexural  mode  and  their 
displacement  is  ideally  zero.  The  first  flexural  mode  shape  is  also  symmetric  about  the 
beam  mid-span,  which  allows  further  simplification.  In  this  reduced  form,  only  the 
support  beam  torsional  stiffness,  Kj,  contributes. 

The  natural  frequency  of  the  first  flexural  mode  of  the  xylophone  bar  is  then 
determined  from  the  lowest  value  of  a  that  causes  the  determinant  to  vanish: 


=  0  (?) 


[sin(a)  +  sinh(a)]  [cos(a)  +  cosh(a)]  0 
[cosh(a)  +  cos(a)]  [-  sin(a)  +  sinh(a)]  - 1 

[sin(a)  -  sinh(a)]a  [cos(a)  -  cosh(a)]a  -  K 

0  0  cos(ra) 

0  0  -  cos(ra) 


0 

0 

-2a 


0 

-1 

-K 


-  [sin(ra)  +  sinh(ra)]  cosh(ra) 
[sin(ra)  -  sinh(ra)]  cosh(ra)J| 


The  parameters  K  and  r  in  Eq.  (7)  are  given  by 

K  =  ^rk  i  =0.2241, and  r  =  —  -1*1.232, 
El  1  2LX 


where  the  length  L\  is  the  distance  from  the  ends  of  the  xylophone  bar  to  the  center  of  the 
supports  and  the  torsional  spring  stiffness,  Kj,  is  derived  assuming  the  xylophone  bar 
applies  a  concentrated  torque  at  the  mid-span  of  the  support  arm.  Hence, 


Kt  = 


4 JSG 

Ls 


(8) 


where  G  =  65  GPa  is  the  shear  modulus  for  polysilicon,  Ls  is  the  length  of  the  support 
arm,  and  Js  is  a  torsional  constant  that  relates  the  angle  of  twist  to  the  applied  torque  and 
is  given  as 


Js=atf 


1  „  „  b 

f  b4 

—  0.210- 

1  < 

3  a 

l  12a4 JJ 

(9) 


The  constants  a  and  b  are  the  width  and  thickness  of  the  support  beam,  respectively.  In 
this  context,  the  thickness  is  defined  as  the  lesser  of  the  two  cross-sectional  dimensions. 

The  smallest  non-zero  a  satisfying  Eq.  (7),  cto,  is  used  to  determine  the  natural 
frequency,  ^Torsion,  of  the  first  bending  mode  including  the  effects  of  the  support  arm 
torsional  resistance: 
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The  form  of  this  expression  is  very  similar  to  that  of  the  free-free  beam,  where  the  first 
bending  mode  frequency,  ^ree  is  given  by 


f Free 


,  with  /?«  4.73 


(11). 


The  calculated  ratio  of  the  these  two  frequencies  is  simply 


/r 
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f Free 
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(12) 


The  measured  natural  frequencies  for  each  of  the  three  variants  of  the  polysilicon 
xylophone  bar  are  plotted  in  figure  6.  The  theoretical  frequencies,  as  outlined  above,  are 
plotted  as  the  continuous  curve.  The  average  percent  error  in  the  theoretical  values  is  less 
than  1%.  The  tendency  to  underestimate  the  actual  resonance  frequencies  may  be  a 
consequence  of:  (a)  assuming  the  applied  torque  to  be  concentrated  at  the  mid-span  rather 
than  distributed  over  the  50  pm  support  arm;  (b)  the  assumptions  made  in  deriving  the 
torsional  spring  stiffness  of  the  support  arms;  or  (c)  variations  in  manufacturing 
tolerances  during  processing. 


Figure  6.  Resonance  frequency  as  a  function  of support  arm  width  for  500x  50  pm 

polysilicon  xylophone  2pm  thick.  The  solid  line  is  the  theoretical  curve.  The  circles 

represent  measured  data. 


A  plot  of  the  output  response  of  the  polysilicon  xylophone  magnetometer  with  4 
pm  support  arm  widths  as  a  function  of  impressed  magnetic  flux  density  up  to  150  pT  is 
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shown  in  figure  7.  The  superimposed  line  is  a  least-squares  fit  to  the  data,  and  it  is 
apparent  that  the  magnetometer  response  is  linear  over  the  range  of  magnetic  flux  density 
shown.  An  estimate  of  the  ultimate  noise  floor  of  the  xylophone  magnetometer,  based  on 
thermomechanical  noise  (Brownian  motion  due  to  the  Langevin  force)  and,  to  a  lesser 
extent,  Johnson  noise  (current  noise  that  couples  with  the  magnetic  field),  yields  a  value 
of  the  order  of  100  pT*AWHz.  This  estimate  calls  out  the  very  important  role  that  the 
amplitude  of  the  drive  current  plays  in  determining  the  ultimate  field  sensitivity  of 
xylophone  bar  magnetometers.  The  decreased  sensitivity  of  these  initial  polysilicon 
devices,  when  compared  to  conventional  (5  x  0.5  x  0.09  mm)  devices  chemically  milled 
from  CuBe  foil,  is  identified  as  being  a  direct  result  of  the  high  sheet  resistivity  of  the 
polysilicon  material  (10  Q/square)  and  its  lower  current  carrying  capability. 


Figure  7.  Output  response  of  a  500x50x2  fjm  xylophone  magnetometer  as  a 
function  of  impressed  magnetic  flux  density 


It  is  well  known  that  air  pressure  greatly  influences  the  performance  of  MEMS 
resonators  and  an  example  of  this  effect  on  the  xylophone  with  4  pm  wide  arms  is  shown 
in  figure  8.  For  this  measurement,  the  magnetic  flux  density  was  on  the  order  of  10  pT 
(Earth  +  other  background  fields)  and  the  xylophone  current  was  22  pA.  The  test 
chamber  pressure  was  varied  between  15  mTorr  and  85  mTorr  and  the  resulting  Q-factors 
were  determined  to  be  21,950  and  2,780,  respectively.  The  9  Hz  decrease  in  resonance 
frequency  at  the  lower  pressure  is  believed  to  be  due  to  thermal  effects  (lower  stiffness). 
Moderate  temperature  sensitivity  of  the  resonance  frequency  has  also  been  observed  at 
constant  pressure  by  varying  the  intensity  of  a  broad  illumination  source  or  by 
repositioning  the  mW  beam  deflection  laser.  There  has  also  been  evidence  that 
mechanical  energy  losses  at  resonance  contribute  enough  thermal  energy  to  measurably 
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lower  the  resonance  frequency  of  a  device. 


78.00  78.05  78.10  78.15  78.20  78.25  78.30 

Frequency  (kHz) 


Figure  8.  Response  of 500x50  pm  xylophone  with  4  pan  supports  at  15mTorr 

(squares)  and  85  mTorr  (circles). 


Figure  9.  SEM  micrographs  of  MUMPs30  poly2  devices  with  gold  deposited  on 

the  electrodes  (a),  as  partial  stripes  on  the  support  arms  (b),  as  full  stripes  on  the 
support  arms  (c),  and,  as  stripes  along  the  xylophone  bar  between  the  support  arms 
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Various  designs  to  decrease  the  resistance  of  polysilicon  devices  were  laid  out  and 
manufactured  on  the  MUMPs30  production  run.  In  these  new  designs  the  drive  current 
to  the  magnetometer  is  fed  from  both  sides,  giving  an  immediate  doubling  of  the  drive 
current  capability.  In  addition,  magnetometers  were  manufactured  from  polyl/poly2 
layers  and  the  use  of  gold  stripes  deposited  on  various  parts  of  the  xylophone  was 
examined.  Also  included  on  the  die  were  various  geometries  of  polyO  and  poly2 
capacitive  plates.  SEM  micrographs  of  released  and  packaged  MUMPs30  poly2  devices 
with  gold  deposited  on  the  electrodes,  as  partial  stripes  on  the  support  arms,  as  full  stripes 
on  the  support  arms,  and,  as  stripes  along  the  xylophone  bar  between  the  support  arms, 
are  reproduced  in  figure  9. 

Characterization  of  these  devices  showed  that  the  polyl/poly2  500  x  50  x  3.5  pm 
bars  with  gold  deposited  on  only  the  electrodes  [figure  9(a)]  have  resonance  frequencies 
of  165  kHz  and  Q-factors  of  between  30,000  and  40,000  at  low  pressures.  The  same  bars 
with  full  gold  strips  on  the  support  arms  [figure  9(c)]  have  increased  resonance 
frequencies  to  168  kHz,  due  to  the  increased  stiffness  of  the  support  arms,  and  Q-factors 
of  10,000.  The  bars  with  the  gold  strips  along  the  xylophone  bar  [figure  9(d)]  have 
decreased  resonance  frequencies  of  approximately  164  kHz,  due  either  to  mass  loading  of 
the  bar  or  to  the  decreased  effective  Young’s  modulus,  and  Q-factors  in  excess  of  10,000. 
From  these  results  and  from  those  from  similar  devices  manufactured  on  subsequent 
MUMPs  runs,  it  is  concluded  that  although  the  use  of  gold  increases  the  current  carrying 
capacity,  there  is  a  concomitant  decrease  in  Q-factor.  As  a  result  there  is  little,  if  any, 
increase  in  overall  sensitivity. 

3.  HYBRID  MAGNETOMETER 

Before  undertaking  detailed  design,  a  study  was  made  of  the  Analog  Devices’ 
BIMOSIIE(iMEMS)  process.  This  process  integrates  an  analog  IC  process  with 
MOSFETs,  BITS,  diodes  and  resistors  with  MEMS  to  create  electro-mechanical  sensors. 
This  process  is  used  to  fabricate  the  ADXL  family  of  accelerometers  sold  by  Analog 
Devices.  The  accelerometers  use  interdigitated  electrodes  to  form  a  capacitive  sensor 
whose  movements  are  restricted  to  a  plane  parallel  to  the  die  surface.  Although  the 
process  is  well  suited  for  this  application,  it  has  been  highly  optimized  and  designs  that 
deviate  considerably  from  this  method  of  sensing  movement  are  not  easily  implemented. 

According  to  the  design  rules  obtained  via  the  MEMS  Exchange,  the  sensor  area 
is  defined  by  the  mask  layer  MOAT.  In  this  region,  all  structural  oxides  will  be  removed 
to  allow  for  release  of  the  mechanical  components  of  the  sensor.  As  a  direct 
consequence,  the  only  elements  that  are  allowed  in  the  sensor  region  are  structural 
polysilicon  for  beams  (BEAM),  polysilicon  anchors  (ANCH),  spacer  oxides  (MOBE)  and 
an  implanted  wiring  layer  (EMIT).  Therefore  the  only  method  of  constructing  the 
xylophone  based  magnetometer  would  be  to  use  BEAM  as  the  vibrating  bar  over  an  area 
of  EMIT  to  form  the  plate  of  a  sense  capacitor. 

This  technique  leads  to  several  problems  for  this  design.  The  capacitive  coupling 
between  the  BEAM  and  plate  implemented  in  the  EMIT  layer  will  not  be  as  strong  as  the 
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present  Cronos  based  design.  The  BIMOSIIE  design  manual  specifies  that  the  structural 
poly  is  1.6pm  above  the  silicon  nitride  layer  that  coats  the  die.  The  nitride  thickness  was 
not  specified.  However,  the  1.6  pm  gap  is  already  considerably  larger  than  the  0.75  pm 
gap  available  between  the  poly2  and  polyl  layers  on  the  Cronos  MUMPs  process. 
Therefore,  the  capacitance  between  the  sense  plate  and  the  beam  will  decrease.  More 
significantly,  the  percent  change  in  the  capacitance  as  the  bar  is  deflected  will  be  reduced, 
assuming  similar  deflections  for  both  devices. 

A  second  problem  arises  by  implementing  sensing  capacitor  plates  in  the  EMIT 
layer.  The  EMIT  layer  is  an  N+  implanted  region  within  the  p-type  silicon  substrate. 
This  forms  a  reverse  biased  diode  between  wiring  drawn  in  EMIT  and  the  substrate. 
Consequently,  there  is  a  large  junction  capacitance  between  the  plates  fabricated  in  the 
EMIT  layer  and  the  silicon  die,  which  is  electrically  connected  to  the  circuit  ground.  If 
the  coupling  area  of  the  structure  is  100x50  pm2,  the  stray  capacitance  to  ground  is  630fF 
—  at  least  20  times  larger  than  the  sensor  capacitance.  This  large  stray  capacitance  would 
significantly  impact  the  first  stage  amplifier  from  the  sensor. 

Finally,  it  is  desirable  for  the  xylophone  to  present  a  low  impedance  to  the  drive 
circuitry.  The  POLY  material  used  to  form  the  bar  is  lOQ/sq,  which  is  comparable  to  the 
present  design.  However,  the  stray  resistance  in  the  wires  and  contacts  to/from  the  bar  is 
quite  large.  The  design  rules  do  not  allow  POLY  or  METAL  wires  to  cross  the  boundary 
between  the  circuit  and  sensor  areas.  Therefore,  we  would  be  forced  to  use  wires  made  in 
the  EMIT  layer  to  interconnect  the  polysilicon  BEAM  and  the  circuitry.  The  design 
manual  specifies  that  the  resistance  of  a  single,  minimum  size  contact  between  the 
polysilicon  BEAM  and  the  EMIT  wire  is  1.2  kO.  Additionally,  the  resistance  of  the 
EMIT  layer  wire  is  14  Q  /sq.  The  total  DC  impedance  will  then  be  a  few  kilo-ohms  -  at 
least  an  order  of  magnitude  larger  than  the  present  design.  Additional  AC  loading  will 
also  be  present  because  of  the  junction  capacitance  of  the  wiring.  Running  the  circuits  at 
a  lower  drive  current  is  not  an  acceptable  means  of  addressing  the  problem  since  the 
lower  current  results  in  less  sensitivity  to  applied  magnetic  field. 

The  Analog  Devices  BIMOSIIE  process  may  seem  attractive  for  use  in  creating 
MEMS  devices  with  integrated  signal  conditioning  and  control  electronics  on  a  single 
die.  However,  the  level  of  optimization  that  has  been  performed  to  construct  the  ADXL 
accelerometers,  does  not  allow  for  simple  migration  of  alternate  sensor  designs.  It  is 
concluded  that  the  poly  silicon  xylophone  bar  magnetometer  project  would  benefit  more 
by  taking  a  hybrid  approach  to  the  system  design.  The  hybrid  would  consist  of  two  dice, 
one  containing  the  mechanical  structures  and  a  second  die  containing  the  associated  drive 
and  sense  electronics.  This  initial  design  of  the  capacitive  sense  using  a  charge  sampling 
circuit  was  fabricated  using  the  AMIS  0.8  pm  process  through  MOSIS. 

The  purpose  of  the  charge  sampling  circuit  is  to  measure  the  deflection  of  the 
magnetometer  transducer  shown  in  figure  10.  A  current,  lac,  is  fed  from  endl  to  end2.  In 
the  presence  of  a  magnetic  field  the  Lorentz  Force  causes  the  beam  to  deflect.  When  the 
beam  bends  down  in  the  center,  the  middle  plate  capacitor  increases  and  the  side  plate 
capacitors  decrease.  Likewise,  when  the  beam  bends  up  the  middle  plate  capacitor 
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decreases  and  the  side  plate  capacitors  decrease.  The  peak  capacitance  change  to  a  lpT 
field  is  approximately  0. 1  %. 


Figure  10.  Magnetometer  transducer 


Figure  1 1  shows  a  circuit  model  for  the  transducer.  Note  that  the  top  plate  of  each 
capacitor  is  also  the  beam.  Resistance  has  been  added  to  model  the  resistance  of  the  beam 
material  as  shown.  Because  the  top  plates  are  connected  to  the  beam,  they  cannot  be  used 
as  the  summing  node  in  a  switched  capacitor  circuit.  Thus,  two  circuits  were  investigated 
that  uses  the  bottom  plates  in  the  summing  node.  The  first  is  a  circuit  that  uses  charge 
redistribution  to  determine  the  difference  in  the  middle  and  side  capacitors.  The  other 
more  robust  design  uses  active  common  mode  subtraction  in  determining  the  capacitance 
difference.  Both  circuits  are  described  and  an  explanation  given  why  the  active  common 
mode  subtraction  was  chosen. 


Figure  11.  Transducer  model 


The  charge  redistribution  switched  capacitor  network  is  shown  in  figure  12.  The 
clocks  are  shown  in  figure  13.  Clocks  phij  and  phi2  are  non-overlapping  clocks.  Clocks 
ending  with  a  ‘d’  are  slightly  delayed  versions.  For  example,  phi2d  is  slightly  delayed 
from  phi2.  In  the  case  of  a  delayed  clock,  care  is  taken  to  ensure  the  delayed  clock  does 
not  overlap  the  next  clock  phase.  For  example,  phi2d  does  not  overlap  phil.  In  this  circuit 
phii  is  split  in  to  two  phases  phi  ,.x  and  phi,.Y. 
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Assume  that  the  supply  rails  Vdda  and  vssa  are  at  the  5V  and  OV  respectively.  Also, 
the  common  mode  voltage,  vcm,  is  2.5V.  The  positive  reference  voltage,  vrefP,  is  4V.  The 
negative  reference  voltage,  vrefm,  is  1  V.  The  current  source  Iac  can  sink  and  source 
current  from  vrefP,  vcm,  and  vrefm. 
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Figure  13.  Clocks  for  charge  redistribution  switched  capacitor  network 


During  phi  j.x  the  middle  plate  capacitor  is  charged  to  a  positive  reference  voltage 
'  “vrefP  -  Vcm”*  During  phii.y  the  side  plate  capacitors  are  charged  to  a  negative  reference 
voltage  “Vrefm  -  vcm”*  For  simplicity,  ignore  the  voltage  across  the  beam,  “Iac  *  Rbeam”. 
More  on  the  effect  of  this  voltage  drop  will  be  explained  in  detail  below. 

During  phi2  the  top  plates  are  connected  to  vcm  and  both  the  bottom  plates  are 
connected  to  the  summing  node  of  the  inverting  amplifier.  The  charge  on  the  middle  and 
side  capacitors  redistribute  and  cancel  each  other  out.  If  the  middle  plate  has  more 
charge,  then  the  virtual  ground  at  the  summing  node  forces  this  charge  to  flow  across  the 
integrating  capacitor,  Qp,  and  increase  the  voltage  of  node  out.  Likewise  if  the  middle 
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plate  has  less  charge  than  the  side  capacitors,  then  charge  is  removed  from  Cjp  and 
voltage  of  node  out  decreases.  The  charge  on  Qp  is  never  cleared.  Thus,  this  circuit  is  a 
discrete  integrator  with  equation: 


where 


or 


out(n)  =  out(n  - 1)  +  ■ 


rjf  _  v  _  v 

*  Vrejp  vcm 


out(n )  =  out(n  —  1)  + 


— (Vcm  ~  Vrefm  ) 

AC(n  - 1)  Vref 


C.„ 


2  ’ 


Previously  the  voltage  drop  across  the  beam  “Iac  *  Rbeam”  was  ignored.  Analysis 
of  the  circuit  reveals  that  this  voltage  will  be  sampled.  Fortunately,  the  voltage  will  not  be 
integrated.  Because  the  voltage  across  the  beam  is  present  with  the  same  polarity  during 
both  phii  and  phi2,  the  circuit  can  be  modeled  as  a  continuous  time  amplifier  as  shown  in 
figure  14.  The  equation  for  this  voltage  is: 

out(n)  =  jf-RbeamxIac(n) 

Thus,  the  beam  voltage  will  enter  into  the  measurement.  However,  it  is  not 
integrated  like  the  signal  AC(n).  This  along  with  the  fact  that  Rbeam  is  approximately  20 
ohms  should  help  the  beam  voltage  from  being  a  problem  in  the  sensing  circuit. 


The  first  problem  with  the  proposed  switched  capacitor  circuit  above  is 
converting  it  to  a  fully  differential  circuit.  In  order  to  convert  it  the  transducer  must  be 
modified  to  have  twice  as  many  capacitors.  A  fully  differential  design  is  required  in 
order  to  minimize  the  effects  of  charge  injection  and  common  mode  interference. 
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Another  problem  with  this  design  is  that  each  transducer  capacitor  is  floating 
during  half  of  phij.  Ideally  when  the  capacitor  is  not  connected  it  will  store  its  previous 
voltage.  However,  due  to  stray  capacitances  this  will  not  be  the  case. 


Figure  15  shows  the  network  during  phii.Y.  During  this  phase  the  side  plates  are 
charged  to  the  negative  reference  voltage.  A  parasitic  capacitance  of  lOfF  is  shown  to 
illustrate  the  problem  with  the  floating  middle  plate  capacitor.  The  parasitic  capacitor  will 
cause  the  middle  plate  capacitor  to  sample  the  vrefm  voltage.  While  it  may  only  sample 
1/1 0th  of  the  voltage,  this  is  significant  since  the  capacitance  change  being  measured  is 
only  1/1000*  of  the  resting  capacitance. 

A  more  robust  capacitive  sensing  network  is  shown  in  figure  16.  The  clock 
scheme  is  shown  in  figure  17.  Unlike  the  circuit  in  figure  12,  this  circuit  is  fully 
differential  without  adding  more  capacitors  to  the  transducer.  This  circuit  uses  an 
amplifier  to  remove  the  common  mode  charge  rather  than  charge  redistribution.  This  is 
done  by  the  use  of  the  amplifier  with  two  inverting  inputs.  Details  of  this  amplifier  are 
shown  in  figure  18. 
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During  phij  the  amplifier  is  configured  as  a  unity  gain  buffer  and  the  middle  and 
side  plates  of  the  transducer  are  held  at  vcm.  During  phi2  the  amplifier  continues  to  hold 
the  plates  at  vcm  by  subtracting  of  the  common  mode  charge  on  the  transducer  capacitors 
through  the  capacitors  Cvcmp  and  CVCmm-  The  fully  differential  amplifier  (figure  19)  moves 
the  difference  charge  to  the  integrating  capacitors  Cjp  and  Cim.  Thus,  the  discrete  time 
equation  of  the  output  voltage  is: 

out(n)  =  outp(n )  -  outm(ri)  =  out{n  - 1)  +  Vref  — ,  C,  -  Cip  =  Cim 

The  differential  amplifier’s  output  common  mode  voltage  is  held  at  vcm  by  the 
switched  capacitor  common  mode  feedback  circuit  in  figure  20.  Both  amplifiers  have  a 
nominal  bias  current,  ibias,  of  5pA. 


Figure  1 7.  Clocks  for  active  common  mode  subtraction  network 


The  active  common  mode  subtraction  network  has  several  advantages  over  the 
charge  redistribution  network.  First,  it  is  fully  differential.  Thus,  it  will  have  better 
rejection  of  charge  injection  and  interference.  Second,  it  has  a  simpler  clocking  scheme 
and  does  not  have  the  floating  capacitance  problem. 


This  circuit  was  verified  with  SPICE  simulations  and  was  fabricated  using  the 
AMIS  0.8  pm  process  through  MOSIS. 
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The  test  setup  used  to  evaluate  the  CMOS  capacitance  sensing  circuit  is  shown  in 
figure  21.  The  capacitive  sensing  integrated  circuit  is  mounted  on  test_board3.  The 
VARTS  I2C  interface  board  and  VARTS  board  contain  programmable  logic  and 
interface  logic  to  allow  data  collection  and  communication  between  the  PC  and 
test_board3.  An  external  power  supply  is  used  to  provide  supply  voltages  for  test_board3 
and  other  electronics  external  to  the  PC. 


Two  digital  clock  signals,  ACLK  and  COMP_CLK,  are  used  to  control  the 
capacitive  sensing  circuit  timing.  A  pseudo-random  digital  signal  (DITHER)  is  applied 
to  the  circuit  to  help  reduce  in-band  tones.  The  value  of  the  dither  signal  changes  for 
each  cycle  of  the  COMP_CLK  clock  signal.  The  sensed  capacitance  undergoes  a  1-bit 
analog-to-digital  conversion  and  the  digital  output  is  updated  once  during  each  cycle  of 
the  ACLK  clock  input  signal.  A  timing  diagram  of  the  capacitive  sensing  circuit  digital 
interface  is  shown  in  figure  22. 


Since  the  interface  between  to  the  CMOS  electronics  is  programmable,  two  test 
cases  were  implemented. 
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In  case  (1),  the  phase  difference  (d®)  between  the  comparison  clock  test  signal 
(COMP_CLK),  as  it  leaves  the  VARTS  Board,  and  the  switched-capacitor  (S-C)  circuit 
clock  test  signal  (ACLK)  is  approximately  90°.  In  case  (2),  d®  is  approximately  180°. 

The  data  is  retrieved  and  stored  in  a  PC  using  a  Matlab  program  which  controls 
the  I2C  card.  Once  retrieved,  the  frequency  content  of  the  data  can  be  plotted  and 
additional  post-processing  and  filtering  implemented.  Post-processing  filters  are  used  to 
remove  the  high-frequency  quantization  noise  above  the  frequency  of  interest. 

In  case  (1),  the  signal-to-noise  ratio  (SNR)  was  observed  to  be  approximately  6- 
7dB.  In  case  (2),  the  signal-to-noise  ratio  (SNR)  was  observed  to  be  approximately  8- 
9dB.  Test  conditions  for  each  case  are  shown  in  Table  1 . 


Test  Case 

1 

MmSHm 

K§M§jag| 

42.46189 

42.46 

2.8 

Case  (1)  -b 

42.45671 

42.45-42.46 

2.8 

Case  (1)  -  c  (No  laser) 

- 

- 

- 

Case  (2)  -  a 

42.45546 

42.45 

2.9 

Case  (2)  -  b  (No  laser) 

- 

- 

Table  1.  Test  conditions 

The  frequency  of  the  xylophone  beam  deflection  was  controlled  by  thermally 
exciting  the  beam  with  a  laser  source  modulated  by  a  Stanford  Research  Systems  Model 
DS345  function  generator.  The  beam  deflection  was  measured  using  a  position  sensitive 
detector  (PSD)  and  an  EG&G  7265  DSP  lock-in  amplifier.  There  are  slight  variations  in 
resonant  frequency  for  case  (1)  and  case  (2).  The  reason  for  this  is  not  fully  understood, 
although  it  is  possible  that  the  DC-field  created  between  the  xylophone  beam  and  the 
capacitive  sensing  plates  is  different  in  each  case. 

The  frequency  content  of  the  unfiltered  output  signal  is  shown  in  figure  23  and 
figure  24  for  both  cases.  Figure  25  and  figure  26  show  similar  data  for  post-filtered 
output  signals. 

The  cause  of  slight  variations  in  resonant  frequency  for  case  (1)  and  case  (2)  is 
not  frilly  understood  from  an  electrical  standpoint,  although  it  is  possible  that  the  DC 
electric  field  created  between  the  xylophone  beam  and  the  capacitive  sensing  plates  is 
slightly  different  in  each  case,  resulting  in  a  change  in  electrostatic  attractive  forces  (and 
hence  resonant  frequency).  Alternatively,  this  variation  could  be  due  to  slight  variations 
in  temperature  in  each  test  run. 
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Figure  23.  Unfiltered  output  power  spectrum  of  hybrid  MUMPS  xylophone 

magnetometer  and  CMOS  capacitive  sense  circuit.  Comparison  of  spectra  when  test 
signal  COMPCLK phase  leads  test  signal  ACLK  by  90  ° (Left)  and  when 
COMPCLK  and  ACLK  are  180°  out-of-phase  (Right). 


Figure  24.  Unfiltered  output  power  spectrum  of  hybrid  MUMPS  xylophone 

magnetometer  and  CMOS  capacitive  sense  circuit  in  the  absence  of  xylophone 
deflection.  Comparison  of  spectra  when  test  signal  COMPJCLK  phase  leads  test  signal 
ACLK  by  90°  (Left)  and  when  COMPJCLK  and  ACLK  are  180°  out-of-phase  (Right). 
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4. 


CMOS  MAGNETOMETER 


As  noted  in  section  2  above,  the  sensitivity  of  polysilicon  xylophone  bar 
magnetometers  is  limited  by  the  current  carrying  capacity.  Recent  advances  in 
fabricating  integrated  MEMS  structures  using  post-processing  of  conventional  MOSIS- 
type  CMOS  foundry  structures  [9]  offer  an  attractive  solution  to  increasing  the  current 
carrying  capacity  of  MEMS-based  magnetometers.  CMOS  MEMS  has  been  around  a 
long  time  for  the  fabrication  of  simple  structures. [10]  The  advent  of  deep  reactive  ion 
etching  (DRIE)  has  allowed  the  development  of  complicated  and  high-definition 
structures  with  performances  comparable  with  their  polysilicon  or  silicon-on-insulator 
equivalents.  These  structures  are  composed  of  combinations  of  the  polysilicon  gate 
layer,  the  interconnect  metallizations,  and  isolation  dielectrics  used  in  standard  CMOS 
processing.  Typically,  the  top  metallization  layer  is  used  to  protect  the  electronic 
circuitry  from  the  post-processing  etching  used  to  release  the  micromechanical  structures. 
The  metal  layers  are  electrically  connected  together  in  conventional  comb  finger 
capacitive  drive  and  sense  structures  so  that  the  comb  fingers  function  in  the  same  way  as 
their  homogeneous  polysilicon  or  silicon  counterparts.  However,  metal  layers  in  multi¬ 
conductor  comb  fingers  can  also  be  selectively  interconnected  to  form  two  sidewall 
capacitors  between  comb  fingers  whose  capacitances  change  with  displacements 
perpendicular  to  the  substrate  surface.  Thus,  using  CMOS  MEMS  enables  the 
development  of  complex  micromechanical  structures,  such  as  the  combination  of 
vibrating  bar  xylophone  magnetometer  and  x-,  y-,  and  z-axis  accelerometers  and  angular 
rate  sensors,  with  state-of-the-art  capacitive  sense  performance  on  a  single  chip.  Custom 
designs  can  be  obtained  either  through  the  DARPA-sponsored  Application-Specific 
MEMS  Process  Service  (ASIMPS)  maintained  by  Carnegie  Mellon  University  or  through 
standard  MOSIS-supported  CMOS  foundries. 

A  study  was  undertaken  on  the  expected  performance  benefits  of  MEMS  based 
magnetometers  fabricated  in  CMOS  technology.  Since  a  similar  layer  scheme  is 
available  in  silicon-on-sapphire  (SOS)  technology,  the  potential  benefits  of  alternative 
optical  sensing  was  included  in  the  study.  This  was  considered  relevant  since  capacitive 
sensing  has  limits  and  there  is  a  need  to  develop  alternative  transduction  schemes  that  are 
compatible  with  both  CMOS  and  MEMS.  The  use  of  a  sapphire  substrate  would  allow 
the  development  of  more  sensitive  optical  transduction  schemes. 

The  model  magnetometer  structure  is  depicted  in  plan  view  in  figure  27,  together 
with  an  exaggerated  view  of  the  movement  of  the  bar  relative  to  the  metal  1  and  metal  2 
interdigitated  capacitor  fingers,  assuming  that  there  is  no  residual  stress  in  the  structure 
and  that  in  the  absence  of  movement  the  static  and  dynamic  fingers  are  parallel. 

The  main  bar  is  400  x  19  pm  in  dimension  and  the  supports  are  30  x  3  pm.  The 
interdigitated  capacitors  are  12  pm  long  with  width  and  gap  spacings  of  1.4  pm.  The 
overlap  between  the  static  and  dynamic  is  assumed  to  be  10.5  pm.  There  are  15  sets  of 
interdigitated  fingers  on  the  bar  between  the  ends  and  the  supports  and  37  sets  between 
the  supports. 
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Metal  2 


Figure  2 1.  Plan  view  of  CMOS-MEMS  xylophone  bar  and  exaggerated  side  view  of 

the  movement  of  the  bar  relative  to  the  static  metall  and  metal2  interdigitated  capacitor 

plates 

One  of  the  unknowns  is  the  form  of  the  xylophone  bar  while  vibrating.  Three 
possible  examples  are  shown  in  figure  28.  These  are  quadratic,  the  classical  Timoshenko 
form  of  a  simply  supported  beam  with  a  concentrated  force  applied  to  the  center  of  the 
bar,  and  the  Timoshenko  form  in  the  center  of  the  bar  but  with  straight  sections  tangential 
to  the  bar  from  the  supports  to  the  ends.  This  latter  form  is  possible  from  the  few  laser 
vibrometer  results  we  have  seen  from  ARL,  Adelphi’s  PZT  samples.  The  Timoshenko 
form  with  a  concentrated  force  applied  to  the  center  is  preferred  to  that  of  a  uniformly 
distributed  load  because  the  boundary  conditions  used  do  not  apply  to  the  xylophone 
magnetometer  with  finite  semi-rigid  supports  not  at  the  ends.  It  can  be  seen  that  the 
relative  deflection  between  the  supports  is  quite  insensitive  to  form  but  there  are 
significant  differences  at  the  ends,  affecting  calculations  on  capacitive  changes  with 
deflection. 


a 


Metal  1 


Figure  28.  Possible  profiles  of  the  forms  of  a  vibrating  xylophone  magnetometer 


The  static  deflection  (ds)  of  the  xylophone  bar  with  a  current  I  (Amperes,  A) 
passing  through  it  in  a  magnetic  field  of  B  (Tesla,  T)  is  given  by[5]: 


d.  = 


5B1L\ 
384 EL 


(13) 


where: 


by 

where: 


Lo  =  length  of  bar  between  supports 

E  -  Young’s  modulus 

Im  =  area  moment  of  inertia 

atfl  12  for  rectangular  bar  of  width  a  and  thickness  b 

At  resonance,  and  in  the  absence  of  non-linear  effects,  the  deflection  (dr)  is  given 

dr  =  Qds . (14) 

Q  =  mechanical  quality  factor. 


The  thickness  of  the  individual  layers  (in  pm)  in  typical  CMOS  and  SOS 
processes,  total  thickness,  and  equivalent  Young’s  moduli  (N.m2)  and  density  (kg.m'3) 


CMOS 

SOS 

oxidel  (pm) 

0.9 

0.5 

metall  (pm) 

0.8 

1.0 

oxide2  (pm) 

0.5 

0.5 

metal2  (pm) 

0.8 

1.0 

oxide3  (pm) 

0.4 

0.5 

metal3  (pm) 

1.0 

3.0 

Total  Thickness  (pm) 

4.4 

6.5 

Young's  modulus  (N.m'2) 

4.34E+10 

5.44E+10 

Density  (kg.m-3) 

2.58E+03 

2.63E+03 

Table  2.  Thickness  of  individual  oxide  and  metal  layers,  total  thickness,  and 

effective  Young ’s  moduli  and  densities  for  CMOS  and  SOS  processing 

are  given  in  table  2. 


The  capacitance  as  a  function  of  deflection  for  the  structure  shown  in  figure  27 
and  the  classical  Timoshenko  form  shown  in  figure  28  is  shown  in  figure  29  for  a  metal  1 
thickness  of  0.8  pm.  This  figure  gives  the  capacitance  of  the  sum  of  the  two  sets  of  end 
interdigitated  fingers,  the  sum  of  the  central  set  of  interdigitated  fingers,  and  the  sum  of 
the  two.  The  zero  deflection  values  are  3.8, 4.1,  and  7.9  fF,  respectively.  The  values  for  a 
metall  thickness  of  1.0  pm  are  4.8,  5.2,  and  10.0  fF,  respectively.  The  difference  in 
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capacitance  of  the  total  capacitance  as  a  function  of  deflection  for  the  thickness  of  both 
metall ’s  is  shown  in  figure  30. 
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Figure  29.  Capacitance  of  sum  of  two  sets  of  end  inter  digitated  fingers,  the  sum  of  the 
central  set  of  inter  digitated  fingers,  and  the  sum  of  the  two  as  functions  of  deflection  of  the 

center  of  the  xylophone  bar. 


Figure  30.  Change  in  capacitance  of  sum  of  interdigitated  fingers  as  function  oj 
deflection  of  the  center  of  the  xylophone  bar. 
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The  maximum  current  a  conductor  can  carry  is  limited  by  Joule  heating.  For 
simplicity  the  supports  are  modeled  as  freely  suspended  filaments  in  vacuum.  In  this 
case  the  voltage-temperature  theorem  [11]  states  that: 

v1  =  8f«r,(r)A(r)rfr  =  8  .  .  (15) 

where: 

V  =  applied  voltage 

Ke  =  equivalent  thermal  conductivity 

po  =  resistivity 

Ts  =  temperature  of  support 

T  =  temperature  of  center  of  filament 

Thus  (V  =  IR  and  R  =  poL/ab): 

,  =  l  ffrfr-r,) 

ab\  p0 
or, 

(r_,)  . 

The  absolute  maximum  temperature  corresponds  to  the  melting  point  (Tm)  of  the 
material.  From  past  experience  characterizing  the  polysilicon  devices,  with  support  arms 
typically  45  pm  x  6  (am  x  1.5  pm,  it  has  been  observed  that  the  resonance  curves  degrade 
at  drive  currents  of  between  1  and  5  mA.  The  “fusing”  current  of  such  support  arms  is  43 
mA  for  a  Tm  =  1410  K.  This  is  to  be  compared  with  the  observed  catastrophic  breakdown 
current  of  approximately  25  mA.  The  “effective”  temperature  differential  corresponding 
to  a  5  mA  drive  current  from  equation  (5)  is  15  K.  For  the  purpose  of  the  present  study  it 
is  assumed  that  the  maximum  operating  drive  current  of  CMOS  and  SOS  devices  will  be 
set  by  a  similar  “effective”  temperature  differential.  It  is  further  assumed  that  oxide3  is  a 
perfect  thermal  insulator  such  that  the  temperature  of  metaB  (through  which  the  drive 
current  is  passed)  is  independent  of  the  underlying  layers.  With  these  assumption  the 
maximum  operating  current  for  the  CMOS  device  is  85  mA  and  for  the  SOS  device  is 
250  mA. 

With  these  values  in  place  the  deflection  of  a  CMOS  and  SOS  xylophone 
magnetometer  with  a  mechanical  Q’s  of  10,000  as  a  function  of  magnetic  field  is  shown 
in  figure  31.  From  this  figure  it  can  be  seen  that  the  sensitivities  are  4.6x1 0'3  m.T'1  for 
the  CMOS  device  and  3.4xl0'3  m.T'1  for  the  SOS  device.  For  comparison,  figure  31  also 
gives  the  deflection  of  a  1000x100x2  pm  polysilicon  device  with  Q  =  10,000  and  a  drive 
current  of  100  pA. 


(16) 

(17) 
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Figure  31.  Deflection  as  function  of  magnetic  field for  a  CMOS  xylophone  with 

Q=10,000  and  a  drive  current  of  85  mA  and  a  SOS  xylophone  with  Q=10,000  and  a  drive 
current  of 250  mA.  Also  shown  is  a  1000x100x2  pm  polysilicon  device  with  q=10,000  and 

a  drive  current  of  0.1mA. 


The  signal  to  noise  ratio  in  the  capacitive  sampling  circuit  due  to  capacitor  noise 
charge  is  given  by  [12]: 


(18) 


where: 

S  =  sensitivity  (m.T'1) 

do  =  static  gap 

Co  =  static  capacitance 

OSR  =  over-sampling  ratio. 

Using  the  sensitivities  and  zero  deflection  capacitances  derived  above,  do  =  1.4 
pm,  an  OSR  of  128,  Vref=5  V,  and  kT=  3.6  x  10-21  V2F,  the  amplitude  of  the  magnetic 
field  for  a  6dB  SNR  is: 

\B±\  *  21  nT 

The  change  in  capacitance  of  a  xylophone  magnetometer  as  a  function  of 
magnetic  field  has  been  shown  [12]  to  be: 
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8C 

8B± 


Co 


This  corresponds  to: 


and, 


AC 
A B± 


-2.6x10  5  fF.nT'1  for  CMOS  devices 


AC 
A B± 


-2.4xl0~5  fF.nT-1  for  SOS  devices. 


(19) 


For  CMOS  devices,  a  field  of  21  nT  corresponds  (using  AQ  =  VAC)  to  a 
minimum  detectable  charge  of  2.7xl0-18  Coulombs,  or  16  electrons  (1  electron  = 
1.6x1  O'1 9  C).  Furthermore,  this  method  implies  that  the  theoretical  minimum  detectable 
magnetic  field  is  about  1.2nT  (field  due  to  charge  from  a  single  electron). 


Inspection  of  (6)  suggests  that  doubling  OSR  decreases  the  minimum  detectable 
field  by  -30%.  Increasing  OSR  by  a  factor  of  16,  for  example,  reduces  the  minimum 
detectable  field  to  5.1  nT  and  increases  the  sampling  frequency  to  96  MHz.  Increasing 
the  sampling  rate,  however,  affects  the  noise  analysis  if  the  bandwidth  of  the  amplifier 
used  in  [12]  is  less  than  the  bandwidth  of  the  sampling  network. 

In  practice,  the  minimum  detectable  field  will  also  be  limited  by  the  ratio  of  the 
dynamic  range  of  the  analog  output  to  the  analog  voltage  that  corresponds  to  the 
minimum  detectable  field,  and  by  the  resolution  of  the  device  used  to  measure  the  analog 
voltage.  It  should  be  noted  that  a  sensor  with  a  digital  output  does  not  have  this 
limitation. 


In  general,  the  feasibility  of  increasing  the  over-sampling  ratio  should  be 
investigated.  In  addition,  a  more  complete  analysis  should  include  noise  contributions 
other  than  thermal  noise. 

Among  the  many  optical  detection  schemes  have  been  implemented  to  detect  the 
displacement  of  MEMS  components,  the  most  common  are  the  optical  lever  and 
interferometry.  The  optical  lever  [13]  is  simpler  to  implement  since  it  does  not  require 
the  positioning  of  the  input  and  output  beams  normal  to  the  reflecting  surface.  However 

the  resolution  is  typically  limited  to  approximately  10"2  nm/  VHz .  Fiber  optic 
interferometers  [14],  based  on  the  interference  of  light  between  the  cleaved  end  of  a  fiber 

and  the  MEMS  device,  generally  have  resolutions  of  the  order  of  10-3  nm/  VHz  .  State- 
of-the-art  optical  interferometers  using  transparent  substrates  have  demonstrated 
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sensitivities  of  6xl0'5  nm/ VHz  with  microcavities  [15]  and  2xl0'5 
diffraction  elements  [16]. 


nm/VHz  with 


It  is  expected  that  the  sensitivity  of  an  easily  implemented  interferometric  scheme 
will  be  of  the  order  of  1CT4  nm/VHz.  With  a  detection  bandwidth  of  ,1  Hz  this 
corresponds  to  a  minimum  displacement  of  10'4  nm.  With  the  sensitivities  derived  in 
Section  3,  the  resultant  minimum  magnetic  field  expected  to  be  detected  with  a  400  x  19 
pm  xylophone  bar  magnetometer  manufactured  in  either  CMOS  or  SOS  technology  is 
approximately  0.025  nT. 

Using  a  detection  bandwidth  of  1  kHz,  as  applicable  in  rotating  platforms,  the 
corresponding  minimum  displacement  is  of  the  order  of  3x1  O'3  nm  and  the  minimum 
detectable  magnetic  field  is  approximately  0.8  nT. 

It  is  anticipated  that  a  working  optical  interferometric  detection  scheme  will  be 
easier  to  implement  using,  at  least  in  the  first  instance,  a  hybrid  approach  and  the 
transparent  substrate  of  SOS  technology.  The  transparent  substrate  will  be  a  pre-requisite 
for  a  truly  integrated  solution. 

In  parallel  with  this  study  the  DARPA-sponsored  Application-Specific  Integrated 
MEMS  Process  Service  (ASIMPS)  Alpha-users  course  was  attended  and  a  CMOS-based 
magnetometer  was  designed.  A  photograph  of  a  released  structure,  as  delivered  from 
CMU,  is  shown  in  figure  32.  At  the  time  of  writing,  it  had  not  been  possible  to  obtain 
any  vibrational  characteristics  from  the  structure. 


Figure  32  SEM photograph  of  released  CMOS-MEMS  xylophone  bar 

magnetometers 


5.  ARL(Aberdeen) 


Some  polysilicon  MUMPs  device  structures,  including  xylophones,  cantilevers, 
and  diaphragms,  were  selected  for  shock  testing,  diced  into  5x5  mm  chips,  released, 
and  mounted  and  bonded  into  40-pin  dual-in-line  ceramic  packages.  An  empty  such 
package  was  fitted  with  a  sapphire  window  and  subjected  to  shock  testing  to  verify 
structural  survivability.  The  package  was  placed  in  an  aluminum  housing  and 
encapsulated  in  a  clear  potting  material.  A  “window”  was  fabricated  in  the  top  surface 
by  milling  out  a  small  patch  and  refilling  by  hand  pouring  a  thin  layer  of  epoxy.  This 
technique  allows  for  the  optical  inspection  and  characterization  of  the  sensor  die 
without  having  to  release  them  from  the  epoxy.  The  packages  were  shock  tested  from 
1,000  to  40,000  g's  ,  0.5  to  0.1ms  duration  and  survived  all  shocks  in  this  range. 
Therefore,  the  device  packages  were  fitted  with  sapphire  windows  were  shock  tested  in 
two  orientations  up  to  30  Kg's.  All  of  the  magnetometer  elements  survived  the  shock 
testing.  However,  it  was  observed  that  a  few  of  the  magnetometer  elements  became 
stuck  to  the  capacitor  plates  as  a  result  of  the  shock  testing.  It  was  found  that  the  stuck 
elements  could  be  released  by  subjecting  them  to  a  large  Lorentz  force  with  a  strong 
magnet. 

The  anticipated  increase  in  shock  testing  did  not  materialize  because  of  problems 
in  developing  reliable  eutectic  die  attach  and  vacuum  packaging  technologies.  Work  at 
ARL/WMRD  therefore  concentrated  on  increasing  the  applicability  of  MEMS-based 
magnetometers  in  projectile  applications. 

The  major  objective  of  the  contract  was  to  develop  a  MEMS-based  magnetometer 
for  use  as  an  orientation  (roll  rate)  sensor  on  projectiles.  As  the  program  progressed,  it 
became  obvious  that  the  class  of  sensitive  miniature  magnetometers  being  developed 
would  find  other  important  applications  in  smart  projectiles  and  bullets,  and  a  particular 
utilization  as  a  compass  or  heading  device.  Given  the  recent  Army/DARPA  MOA  for 
Future  Combat  System  (FCS)  where  mature  technologies  must  be  demonstrated  for 
insertion  into  new  weapon  systems,  this  program  could  be  of  direct  application.  The  new 
munition  concepts  will  most  likely  rely  on  higher  performance  MEMS,  and  a  MEMS- 
based  magnetometer  could  be  an  important  part  of  that  solution. 

Historically,  the  use  of  magnetometers  (as  compasses)  to  measure  the  direction  of 
travel  with  respect  to  the  local  Earth’s  magnetic  field  and  to  aid  inertial  measurements 
units  (IMU’s)  for  guidance  and  navigation  have  been  restricted  to  wing-stabilized 
.  platforms.  Since  most  bullets  and/or  projectiles  roll,  even  those  that  are  fin-stabilized, 
classical  strapdown  IMUs,  configured  from  a  constellation  of  accelerometers  and  angular 
rate  sensors,  will  suffer  from  large  inaccuracies  over  long  flight  times  if  not  aided. 
Conventional  angular  rate  sensors  have  extreme  difficulty  with  scale  factor  and  drift 
errors  on  rolling  airframes.  These  errors  are  due  to  attitude  error  buildup  in  roll,  which 
cross-couples  pitch  plane  dynamics  into  the  yaw  plane.  In  addition,  MEMS-based 
angular  rate  sensors  are  at  an  early  stage  of  development  and  can  be  dramatically 
improved  in  performance  if  aided  by  additional  sensors. 
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A  miniature  magnetometer  mounted  on  a  rolling  bullet  would  serve  to 
significantly  increase  the  IMU  accuracy.  A  recent  discovery  at  ARLAVMRD  entitled 
MAGnetic-SONDE  (MAGSONDE)  [17]  -  ARL  patent  pending  -  is  the  theoretical  basis 
for  this  novel  augmentation  of  an  IMU  on  a  spinning  body  containing  rate  sensors  and 
magnetic  sensors.  MAGSONDE  has  been  designed  for  use  in  continuously  rolling 
airframes  for  estimation  of  in-flight  angular  orientation  with  respect  to  the  Earth’s 
magnetic  field.  The  MAGSONDE  is  comprised  of  both  a  magnetic  sensor  and  a 
methodology  that  determines  orientation  from  sensor  phase  measurements.  Importantly, 
sensor  scale  factor  and  gain  variations  will  not  affect  MAGSONDE  performance. 

Although  devices  responsive  to  Earth’s  magnetic  field  have  long  been  used  for 
heading  estimation,  MAGSONDE  is  a  new  and  unique  technology.  It  differs  from  all 
other  known  systems  giving  orientations  with  respect  to  a  magnetic  field  in  that  those 
systems  use  one  or  more  of  four  basic  measurements  to  determine  orientations:  1)  field 
strength  along  a  sensor  axes;  2)  relative  field  strength  along  multiple  sensor  axes;  3)  rate 
of  change  of  field  strength  along  a  sensor  axis;  and,  4)  relative  rates  of  change  along 
sensor  axis.  In  every  case,  analysis  of  the  measurements  is  premised  on  some  evaluation 
of  a  component  of  the  magnetic  field  along  a  sensor  axis  and  requires  prior  knowledge  of 
the  field  and/or  accurate  sensor  gain  calibration.  Making  angular  measurements  using 
MAGSONDE  requires  only  the  magnetic  sensor(s)  to  identify  the  times  at  which  there  is 
no  magnetic  field  along  the  sensor  axis  since  MAGSONDE  determines  orientation  from 
relative  phase  information  in  the  sensor  outputs  and  is  amplitude  independent.  This 
feature  is  important  for  several  reasons:  1)  No  knowledge  of  the  field  strength  is  required, 
2)  manufacturing  tolerances  of  scale  factor  and/or  gain  variations,  which  are  the  bane  of 
many  inertial  sensors,  have  no  effect  on  MAGSONDE  performance,  and  3)  only  scalar 
operations  are  required  for  the  angular  measurements. 

MAGSONDE  could  be  combined  with  additional  field  related  data  sources  to 
yield  :  Projectile  Orientation  in  Navigation  Terms  (POINTER)  [18]  -  ARL  patent  pending 
-  such  that  direct  measurement  of  the  angles  in  a  dual-field  environment  allows  a  derived 
inertial  orientation,  eliminating  integration  of  the  rate  sensor  data  which  is  impracticable 
over  long  periods  of  flight  time  on  spinning  bodies.  Examples  of  additional  fields  are 
sunlight,  gravity,  AGPS  RF  fields,  passive  radiometry,  millimeter  wave  radar,  and 
ground-based  and  flight  body-based  RF  telemetry  fields.  To  make  this  concept  work  the 
MAGSONDE  application  requires  the  development  of  miniature,  low-cost,  bias-free, 
MEMS-based  magnetometers. 

Measurements  of  the  angle  between  the  Earth’s  magnetic  field  and  a  projectile,  as 
illustrated  in  figure  33,  can  be  used  to  track  its  angular  motion.  On  a  spinning  projectile, 
the  orientation  of  the  sensors  is  continuously  changing  with  respect  to  the  field,  unless  the 
spin  axis  is  parallel  to  the  field.  In  this  case  the  measurement  can  be  made  from  phase 
information  in  the  signals  from  appropriately  located  magnetometers.  Roll  rate  with 
respect  to  the  magnetic  field  also  can  be  obtained.  MAGSONDE  consists  of  both  the 
magnetic  sensor  and  the  data  reduction  methodology  that  leads  to  this  measurement 
capability. 
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MAGNETIC  FIELD 


Figure  33.  Illustration  of  angle  (<Jm)  between  Earth ’s  magnetic  field  and  a  projectile 


The  sensor  output  (Os),  illustrated  in  figure  34  for  two  sensors  with  different 
orientations,  depends  on  the  strength  of  the  magnetic  field  (I  Ml ),  the  projectile  angle  in 
that  field  (oM),  the  angle  of  the  sensor  with  respect  to  the  spin  axis  (Xs),  the  roll 
orientation  of  the  projectile  (<J>s)  and  the  sensor  gain  (Gs),  according  to  the  equation: 


Os  =  (Gs)|M|[cos(4)cos(<tm)  +  sin(4)sin(Osin(^)] 


(20) 


Solving  for  <|>s  at  the  zero  crossings  gives  the  projectile  roll  angles  at  which  the 
respective  sensors  are  orthogonal  to  the  magnetic  field.  These  angles  are  a  function  of 
the  sensors’  orientations  on  the  projectile  and  the  projectile’s  orientation  in  the  field,  as 
illustrated  in  figure  35,  and  as  expressed  in  the  following  equation: 


tps  ~ 


sm 


-cos(4)cos  (Ojf) 
sin(4s)sin(o-M) 


(21) 
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are  unequal  and  non-supplementary  is  a  discriminant  of  the  projectiles  orientation  in  the 
field: 


This  is  illustrated  in  figure  36.  The  ambiguity  arising  from  the  symmetry  of  this 
ratio  about  aM  -  0  is  easily  resolved  by  checking  the  parity  of  the  field  along  Si  when  S2 
is  orthogonal  to  the  field. 


The  in-flight  data  will  not  give  sensor  roll  angles  at  zero  crossings  directly  but  the 
times  at  which  these  crossings  occur,  as  illustrated  in  figure  34.  These  times  will  be  used  to 
estimate  the  ratio: 
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(23) 


t(0  -  40 
WJ  -  Wj 


If  for  four  consecutive  zero  crossings  (i.e.  a  roll  cycle),  the  roll  rate  and  projectile 
angle  in  the  field,  ctm,  are  constant,  the  ratio  computed  from  crossing  times  is  equivalent 
to  that  computed  from  roll  angles.  Though  these  constraints  are  seldom  completely  true 
in  practical  flight,  the  estimation  errors  are  usually  small  for  reasonably  stable 
trajectories. 

The  roll  rate  is  estimated  by  numerical  differentiation  of  the  calibrated  crossing 
positions  using  the  measured  crossing  times.  When  ctm  is  near  90°  and/or  the  pitch  and 
yaw  rates  are  low  compared  to  the  body  spin  rate,  the  magnetic  roll  rate  and  body  spin 
rate  are  equivalent. 

If  the  projectile’s  spin  axis  is  colinear  with  the  magnetic  vector  (i.e.  ctm  is  either  0 
or  180°),  the  field  strength  along  the  sensor  axis  is  the  same  throughout  a  roll  cycle  and 
the  projectile’s  magnetic  heading  can  be  determined  with  a  parity  check.  For  other 
projectile  headings  relative  to  the  field,  when  X  ^  90°,  the  ratio  O  can  only  be  computed 
when  (90  -  X)  <  orM  <  (90  +  X).  This  bound  arises  from  the  existence  criterion  for  a 
solution  to  the  equation  of  the  roll  angle  at  which  the  sensor  will  be  orthogonal  to  the 
field  of: 


cos(4)cos(o-m) 

sin(4)sin(<rA/) 


<  1 


(24) 


Though  this  requirement  defines  a  range  of  magnetic  aspect  angles  for  which  the 
ratio  can  be  estimated,  this  limitation  usually  can  be  avoided  by  either  an  appropriate 
choice  of  launcher  aiming  or  the  use  of  multiple  sensors.  When  this  is  not  possible,  it  is 
important  to  note  that  this  limitation  would  probably  only  apply  to  a  portion  of  the 
trajectory  of  many  projectiles. 

Although  the  Earth’s  magnetic  field  varies  both  with  location  and  time,  these 
variations  are  regular  and  known.  Moreover,  the  variations  over  the  length  and  duration  of  a 
projectile  trajectory  are  typically  negligible,  excluding  local  anomalies.  Thus,  given 
knowledge  of  the  flight  location,  the  magnetic  field  near  the  Earth’s  surface  can  be  obtained 
from  geodetic  survey  data,  computer  models,  or  direct  measurement.  Simulated  trajectory 
data,  shown  in  figure  37  for  a  6-inch  diameter,  spin-stabilized  artillery  round  at  the 
Transonic  Range  at  Aberdeen  Proving  Ground,  is  then  used  to  estimate  the  nominal 
anticipated  magnetic  heading  history.  The  corresponding  lower  bound  of  error  in  orientation 
estimate  due  the  violation  of  constraints  outlined  above  is  given  in  figure  38. 
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Figure  38.  Lower  bound  of  error  in  orientation  estimate  due  to  violation  of 


In  summary,  MAGSONDE  determines  the  orientation  from  the  relative  phase 
information  in  the  sensor  outputs  and  is  independent  of  the  magnetic  field  strength. 
This  feature  is  important  for  several  reasons; 

1.  No  knowledge  of  the  field  strength  is  required; 

2.  Manufacturing  tolerances  on  scale  factor  and/or  gain  variations  have  no  effect  on 
MAGSONDE  performance;  and, 

3.  Only  scalar  arithmetic  operations  are  required  for  the  computation  of  angular 
orientations. 

In  order  to  implement  these  concepts  in  a  practicable  manner,  one  requires  a 
sensor  whose  output  is  directly  proportional  to  field  strength  (i.e.  a  DC  sensor),  a  fast 
response,  and  a  wide  dynamic  range.  The  sensor  also  has  to  be  small  and  rugged 
enough  to  be  mounted  on  a  projectile.  This  suggests  the  use  of  MEMS-based 
magnetometers  in  general  and  the  JHU/APL  xylophone  magnetometer  in  particular. 
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6. 


ARL(Adelphi) 


At  the  start  of  the  contract  more  than  50  batches  of  the  PZT  solutions  had  been 
prepared  and  spun  coated  on  65  substrates.  PZT  films  (0.25  to  2.0  pm  thick)  were 
crystallized  using  the  ARL  developed  rapid  thermal  annealing  procedure.  Remnant 
polarization  Pr  (pC/cm2)  and  coercive  field  Ec  (kV/cm)  were  measured  to  verify  the  re¬ 
producibility  of  the  piezoelectric  characteristics  of  the  films.  The  average  measured  Pr 
and  Ec  values  were  20  pC/cm2  and  72  kV/cm,  respectively.  The  thickness  uniformity 
variation  was  found  to  be  ±  2%  across  a  4-inch  substrate. 

A  prototype  magnetometer,  designed  for  process  development  rather  than  operational 
performance,  was  designed  and  is  shown  schematically  in  figure  39. 


\ 

Figure  39.  Schematic  plan  (left)  and  side  (right)  views  of  PZT  device 


The  process  steps  include: 


a. 

b. 

c. 


d. 

e. 

f. 
g- 


sol-gel  deposit  PZT  on  Si/Si02/Pt  substrate; 
sputter  the  top  Pt  electrode; 

photolithography  define  the  xylophone  structure,  ion  mill  the  top 
Pt,  dry  etching  the  PZT,  and  ion  mill  the  bottom  Pt; 
selectively  expose  bottom  Pt  at  contact; 

deposit  PECVD  oxide  and  open  up  contact  vias  for  top  and  bottom 
electrodes; 

deposit  and  lift-off  contact  metallization;  and, 

deep  RIE  etch  substrate  and  oxide  from  beneath  the  xylophone. 


An  SEM  micrograph  of  a  xylophone  structure  that  had  completed  all  the  process 
steps  with  the  exception  of  final  substrate  release,  is  shown  in  figure  40. 
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Progress  on  the  final  release  was  hampered  by  problems  with  the  deep  reactive 
ion  etch  system.  However,  work  continued  to  develop  the  process  to  minimize  the 
residual  stress  in  the  Pt/PZT/Pt  structure  and  thus  allow  the  release  of  an  undistorted 
xylophone  bar.  In  these  studies,  different  layer  sequences  were  deposited  on  3”  silicon 
wafers  and  stress  measurements  were  made  using  a  Tencor  FLX-2908  before  and  after 
adding  each  of  the  layers  as  well  as  during  and  after  annealing.  In  the  case  of  multilayer 
films,  the  total  thickness  of  all  layers  at  the  given  process  step  was  considered.  Negative 
values  for  stress  indicated  a  compressive  stress  and  positive  ones  a  tensile  stress. 
PECVD  oxide  and  nitride  dielectrics  were  deposited  at  250°C  and  rapid  thermal  annealed 
(RTA)  annealed  at  750°C  for  60  sec  in  nitrogen.  Top  and  bottom  metal  electrode  layers 
were  sputter  deposited  and  either  RTA  annealed  at  700°C  for  60  sec  or  furnace  annealed 
to  700°C  (ramped  at  15°C/min)  under  nitrogen.  PZT  films  (250  nm,  500  nm,  1000  nm) 
were  spin  coated  (ref  6)  and  the  residual  stress  was  measured  after  crystallization. 

The  deposition  of  the  bottom  electrode  layer  (20nm  Ta/170  nm  Pt)  onto  either 
dielectric  layer  did  not  significantly  change  the  residual  stress  levels.  However,  after 
RTA  annealing  at  700°C  for  60  seconds,  the  residual  stress  became  tensile  at  levels  of 
1070±70  MPa.  This  is  illustrated  in  figure  41.  This  effect  was  examined  in  more  detail 
by  performing  an  RTA  anneal  at  750°C  after  the  dielectric  deposition  and  before  the 
electrode  deposition.  The  result  of  the  combined  process  is  illustrated  in  figure  42  for 
various  thickness  of  dielectric  layer.  It  can  be  seen  that  the  results  from  the  oxide  layers 
■  are  very  consistent.  The  RTA  annealing  after  oxide  deposition  changes  a  compressive 
stress  of  -430±50  MPa  to  a  tensile  stress  of  +50±50  MPa.  The  stress  then  reverts  to 
compressive  after  deposition  of  the  bottom  electrode,  and  then  becomes  tensile  upon 
RTA  annealing  of  the  stack.  The  residual  stress  decreases  with  increasing  oxide 
thickness  and  varies  from  200  to  500  MPa  -  levels  significantly  lower  than  that  shown  in 
figure  41 .  The  behavior  of  the  nitride-coated  samples  was  much  more  varied  and  resulted 
in  residual  tensile  stresses  in  excess  of  100  MPa.  It  is  to  be  noted  that  the  adhesion  of  the 
subsequently  deposited  PZT  improved  significantly  following  annealing. 
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Figure  41  Residual  stress  after  RTA  annealing  of  combined  as-deposited  bottom 

dielectric  and  electrode  layers 


Figure  42  Residual  stress  after  RTA  annealing  of  separate  as-deposited  bottom 

dielectric  and  electrode  layers 
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The  effect  on  residual  stress  of  the  subsequent  deposition  of  a  500  nm  PZT  layer 
and  a  Pt  top  electrode  layer,  followed  by  an  RTA  anneal  at  300°C  for  60  seconds  in  N2,  is 
illustrated  in  figure  43.  In  each  case,  the  deposition  of  the  PZT  (including  the 
crystallization  process)  reduces  the  residual  stress  to  levels  that  remain  essentially 
unchanged  following  the  top  electrode  deposition  and  RTA  anneal. _ | _ 


Figure  43  Residual  stress  after  deposition  of 500  nm  PZT  layer,  Pt  top 

electrode,  and  RTA  anneal 


The  effect  of  PZT  layer  thickness  was  studied  using  substrates  coated  with  600 
nm  thick  PECVD  oxide  layers  and  the  standard  bottom  electrode.  The  general  behavior 
illustrated  in  figure  43  was  observed  with  the  final  increasing  slightly  with  PZT 
thickness.  For  example,  the  residual  stress  after  recrystallizaton  varied  froml83  MPa  for 
a  250  nm  thick  PZT  layer  to  257  MPa  for  a  1000  nm  thick  PZT  layer.  In  each  case,  the 
residual  stress  reduced  by  approximately  70  MPa  following  the  top  Pt  metallization 
deposition  and  300°C  RTA  anneal.  The  final  residual  stress  level  of  140±25  MPa  tensile 
is  to  be  compared  with  a  residual  stress  level  of  approximately  10  MPa  compressive  for 
polysilicon  layers.  Further  studies  were  undertaken  on  examining  the  development  of 
residual  stress  during  the  deposition  and  RTA  annealing  of  multilayers  of  PZT.  A 
description  of  the  treatments  according  to  the  layers  is  given  in  Table  2.  The  stress  here 
was  based  on  the  entire  thickness  of  the  layer(s).  As  the  layers  accumulated,  the  stress 
values  had  less  variation  and  normalized  to  a  tensile  value  of  175  MPa.  Measured 
negative  values  for  stress  indicate  a  compressive  stress  and  positive  ones  a  tensile  stress. 

The  calculated  residual  stress  of  each  individual  layer  is  given  in  Table  3.  The 
as-deposited  condition  of  the  SiCh,  top  and  bottom  Pt  layers  were  compressive  but 
changed  to  the  tensile  after  annealing.  The  PZT  layers  varied  between  113-148  MPa 
depending  on  thickness. 
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Layer 

Avg. 

Stress 

(Mpa) 

Si02 

-341 

RTA  Si02  (700°C/60  sec  in  N2) 

47 

Ta/Pt  bottom  electrode 

-23 

RTA  Ta/Pt  (700°C/60  sec) 

219 

0.22  pm  PZT 

204 

0.44  pm  PZT 

196 

0.66  pm  PZT 

181 

0.88  pm  PZT 

165 

Pt  top  electrode 

163 

RTA  Pt  (350°C/120  sec) 

175 

St.  Dev. 
(5)± 


Table  2. 


Total  accumulated  successive  stress  in  piezoelectric  multilayer  stack. 


Layer 

Avg. 

Stress 

(Mpa) 

St.  Dev, 
(6)± 

Si02 

-341 

7 

RTA  Si02  (700°C/60  sec  in  N2) 

47 

7 

Ta/Pt  bottom  electrode 

-284 

14 

RTA  Ta/Pt  (700°C/60  sec) 

858 

6 

0.22  pm  PZT 

144 

7 

0.44  pm  PZT 

148 

9 

0.66  pm  PZT 

132 

14 

0.88  pm  PZT 

113 

18 

Pt  top  electrode 

-113 

26 

RTA  Pt  (350°C/120  sec) 

-15 

37 

Table  3. 


Calculated  residual  stress  of  each  individual  layer  in  piezoelectric 
multilayer  stack  during  deposition. 


Individual  layers  were  then  removed  starting  from  the  top  of  the  stack.  The 
measured  stress  was  based  on  the  radius  of  curvature  of  the  wafer  before  and  after  the 
specific  layer  removal.  Removal  of  the  layers  was  accomplished  by  ion  milling  for  the  Pt 
layers,  and  wet  chemical  etching  for  the  PZT  layers.  The  results  are  summarized  in  Table 
4,  together  with  details  of  the  etching  steps.  The  stress  values  of  the  removed  layers 
(Table  4)  were  compared  with  those  of  the  individual  layers  (Table  3).  The  greatest 
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residual  stress  was  found  to  be  in  the  in  the  bottom  electrode  (518  MPa)  after  the  PZT 
was  removed  had  the  largest  difference. 


Etch  method 

Etched  layer 
thickness 
(A°) 

Avg. 

Stress 

(Mpa) 

St.  Dev. 
(5)± 

Ion  milling 

1700  top  Pt  + 

-740  PZT 

-9 

9 

Wet  etch 

3500-6800 

PZT 

183 

19 

Ion  milling 

1900  Ta/Pt  + 
-900  Si02 

518 

31 

Table  4  Stress  due  to  etched  layers. 


Advances  in  process  development  enabled  the  release  of  xylophone  bar  test 
structures  from  the  front  side  of  the  wafer  rather  than  from  the  back  side.  Holes  were 
opened  up  in  the  oxide  surrounding  the  xylophone  structures  and  the  silicon  was  removed 
by  use  of  a  40-pm  deep  “isotropic”  etch  in  the  DRIE  system.  This  resulted  in  fully 
released  and  suspended  PZT  xylophone  bar  magnetometer  structures  as  shown  in  figure 
44.  The  isotropic  release  process  using  a  SF6  plasma  performed  well  and  did  not  leave 
any  residue  on  the  top  metallization.  The  xylophone  bars  have  been  electrically  excited 
in  order  to  prove  piezoelectric  actuation  within  the  structure  and  to  locate  the 
fundamental  resonance  frequencies  of  the  bars. 


Figure  44.  Resultant  PZT  xylophone  bar  actuator  with  the  ends  of  the  bar  bowing 

out  of  plane  due  to  the  large  tensile  stresses  in  the  platinum  metallization. 


Real-time  structural  dynamics  measurements  on  a  nominally  1  x  0.1  mm  device 
were  obtained  using  a  laser  Doppler  vibrometer.  The  xylophone  bar,  operating  in  air,  was 
actuated  with  a  1  volt  input  signal  applied  to  the  top  electrode  which  produced  a 
piezoelectric  response.  Results  of  the  three  lowest  frequency  modes  observed  are  shown 
in  figure  45.  The  resonance  frequencies  of  the  torsion  mode  (figure  45a),  sine  wave 
mode  (figure  45b),  and  the  trampoline  mode  (figure  45c)  were  14.84  kHz,  16.03  kHz, 
and  26.47  kHz,  respectively. 


Figure  45  Laser  vibrometer  scans  showing  torsion  mode(a),  sine  wave  mode  (b), 
and  trampoline  mode  (c)  deflections  of  a  xylophone  bar  test  structure. 


A  part-wafer  (approximately  38  x  17  mm  in  size)  containing  several  operational  devices 
was  used  in  an  attempt  to  characterize  their  magnetometer  performance.  Unfortunately, 
the  devices  on  the  wafer  had  been  aligned  45°  to  the  cleavage  planes  of  the  substrate  and, 
as  a  consequence,  could  not  be  reduced  to  a  manageable  size  without  endangering  the 
integrity  of  the  released  xylophone  bars.  Therefore,  the  whole  piece  was  mounted  on  a 
large  header,  which  in  turn  was  adapted  to  fit  into  the  vacuum  test  fixture. 

Initial  characterization  was  undertaken  by  applying  a  1  or  10  mV  ac  signal  across 
the  top  electrode  and  sweeping  the  frequency  from  1  to  100  kHz.  The  results  for  a 
nominally  1  x  0.1  mm  and  a  nominally  0.5  x  0.05  mm  device  are  shown  in  figures  46(a) 
and  46(b),  respectively.  In  these  figures,  the  red  curves  are  the  magnitude  of  the  beam 
deflection  signals  and  the  blue  curves  are  the  phase  outputs.  Also  shown  are  SEM 
micrographs  of  the  two  types  of  device.  It  can  be  seen  that  the  frequency  spectrum  is 
very  complex,  probably  as  a  result  of  the  non-idealized  tether  structure. 

Actual  dimensions  of  the  structures  were  taken  from  an  automated  micrometer 
microscope  stage.  These  dimensions,  together  with  the  thickness,  density,  and  Young's 
modulus  of  each  layer  in  the  structure  resulted  in  expected  frequencies  of  the  first 
torsional  mode  of  7.34  kHz  and  30.13  kHz  for  the  two  sizes.  This  is  to  be  compared  with 
the  observed  frequencies  of  8.31  and  42.4  kHz,  respectively.  The  differences  in 
frequencies  observed  in  the  laser  vibrometer  [figure  45]  and  in  the  beam  deflection 
[figure  46(a)]  is  attributed  to  the  devices  coming  from  different  wafers,  with  the  former 
having  unknown  thickness  of  active  components.  It  is  also  to  be  noted  that  significant 
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differences  were  observed  among  the  different  devices  on  the  same  wafer,  indicating 
variation  in  material  characteristics  across  the  wafer. 
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Figure  46.  Frequency  scans  and  SEM  micrographs  of  the  two  sizes  of  prototype 
PZT xylophone  structures  -  (a)  nominally  1  x  0.1  mm,  and  (b)  nominally  0.5  x  0.05 

mm. 


The  region  of  the  lowest  frequency  mode  in  all  devices  tested  to  date  contained 
two  peaks.  Figure  47  shows  these  peaks  in  some  detail  in  an  experiment  designed  to 
determine  the  response  of  the  prototype  devices  to  a  magnetic  field.  This  problem  was 
compounded  by  the  fact  that  the  only  header  available  was  deep  and  magnetic  and 
effectively  shielded  the  wafer  from  any  impressed  field.  Figure  47(a)  shows  the 
frequency-scan  across  the  two  peaks  with  a  drive  current  applied  to  the  bottom  electrode 
only.  The  response  could  be  due  to  either  a  Lorentz  force  between  the  drive  current  and 
the  (unknown)  magnetic  field  present,  the  piezoelectric  effect  arising  from  the  difference 
in  field  across  the  PZT  layer  from  the  voltage  drop  along  the  bottom  layer,  or  a  mixture 
of  both.  Figure  47(b)  shows  the  frequency-scan  with  the  drive  current  applied  to  both  top 
and  bottom  electrodes  in  the  same  direction.  This  configuration  should  increase  the 
Lorentz  contribution  and  decrease  the  piezoelectric  contribution.  It  can  be  seen  that  the 
overall  magnitude  of  both  peaks  is  decreased  by  almost  two  orders  of  magnitude, 
implying  that  the  device  is  being  driven  primarily  by  the  piezoelectric  response  of  the 
PZT,  and  that  the  ratio  of  intensities  of  the  two  peaks  is  different.  Figure  47(c)  shows  the 
frequency  scan  with  the  drive  current  applied  to  both  top  and  bottom  electrodes  in  the 
same  direction  and  in  the  presence  of  a  strong  (~  lmT)  external  magnetic  field.  There  is 
no  significant  difference  between  figures  47(b)  and  47(c)  -  demonstrating  the  magnetic 
shielding  of  the  package.  (The  features  in  the  figure  47(c)  scan  around  44kHz  are 
instrument  noise).  Figure  47(d)  shows  the  frequency-scan  with  the  drive  current  applied 
to  both  top  and  bottom  electrodes  in  the  opposite  direction.  In  this  configuration  the 


Lorentz  force  component  should  be  significantly  reduced  (not  eliminated  because  of 
different  resistances  of  the  top  and  bottom  Pt  layers)  and  the  piezoelectric  component 
enhanced.  It  can  be  seen  that  the  low  frequency  peak  is  visibly  reduced  and  the  high 
frequency  peak  is  comparable  to  that  shown  in  figure  47(a). 
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Figure  47.  Frequency  scans  of  a  nominally  0.5  x  0.05  mm  device  under  different 
drive  conditions  -  (a)  with  the  drive  current  applied  to  the  bottom  electrode  only;  (b) 
with  the  drive  current  applied  to  both  electrodes  and  in  the  same  direction;  (c)  as  (b) 
but  in  the  presence  of  a  magnetic  field;  and,  (d)  as  (b)  but  with  the  currents  flowing  in 

opposite  directions. 


Two  issues  arose  out  of  this  fabrication  process  development.  One  was  the  non- 
planar  nature  of  the  resultant  released  structure,  while  the  other  was  the  low  mechanical 
Q  associated  with  the  xylophone  bar.  The  lower  than  expected  Q  value  was  tentatively 
attributed  to  an  excessive  tether  width,  and  design  modifications  have  been  incorporated 
into  the  mask  set  to  be  used  in  the  first  prototype  device  fabrication.  The  bowing  of  the 
released  xylophone  bar  was  found  to  result  from  the  stress  gradient  in  the  multilayer  stack 
actuator  with  specific  attention  on  the  high  tensile  stress  of  the  underlying  bottom  Ta/Pt 
metallization  layer.  The  SiC>2  thin  film  underneath  the  Pt/PZT/Pt  actuator  causes  the 
neutral  axis  of  the  beam  to  be  located  below  the  bottom  Pt  metallization.  With  the  large 
tensile  stress  of  the  Pt  located  above  the  neutral  axis  the  beam  is  directed  upward  away 
from  the  substrate  surface. 

In  response  to  concerns  over  lower  Q  and  sensitivity  values  due  to  the  non-planar 
released  beam,  several  fabrication  schemes  were  initiated  to  fabricate  planar  PZT 
xylophone  magnetometers.  The  most  promising  result  so  far  was  achieved  using  a  high 
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tensile  stress  Si3N4  layer  placed  in  a  SiC>2  sandwich  underneath  the  Pt/PZT/Pt  actuator. 
If  the  nitride  film  is  placed  properly  below  the  neutral  axis  of  the  beam,  it  can  offset  the 
highly  stressed  bottom  metallization,  thus  yielding  a  stress  compensated,  planar  beam. 

The  study  focused  on  the  fabrication  of  center  clamped  resonating  bars  with  a 
~500A  PECVD  Si3N4  film  deposited  at  different  depths  within  a  PECVD  Si02  sandwich. 
Table  5  lists  the  five  different  fabrication  schemes  that  were  investigated  as  well  as  the 
residual  stress  associated  with  each  film  stack  as  measured  after  a  nitrogen  anneal  to 
remove  excess  trapped  hydrogen.  The  resultant  beams  revealed  a  transition  from  bowing 
in  toward  the  substrate  to  bowing  away  from  it. 


Si02 

Si3N4 

sio2 

Residual  Stress 

(A0) 

(A0) 

(A°) 

(Mpa) 

1000 

500 

5500 

105 

2000 

500 

4500 

120 

3250 

500 

3250 

50.8  | 

4500 

500 

2000 

13.5 

5500 

500 

1000 

-8.5 

Table  5.  Dielectric  structure  of  center-clamped  PZT  resonating  bars  and  the 

residual  stress  in  the  dielectric  stack  after  a  700  °C  anneal  in  N2. 


A  mask  set  for  the  initial  prototype  PZT  xylophone  bar  magnetometers, 
incorporating  modifications  identified  during  the  initial  development  phase,  was 
designed.  This  includes  design  changes  to  the  tether  width  so  that  larger  Q  values  could 
be  realized.  In  addition,  the  mask  set  was  meticulously  designed  to  ensure  that  the  tether 
length  could  be  maintained  with  specific  dimensions  after  release.  Another  feature  that 
became  apparent  during  initial  development  was  the  difficulty  of  dicing  the  wafer  without 
damaging  the  released  structures.  The  dicing  was  addressed  by  designing  a  deep  reactive 
ion  etched  (DRIE)  trench  to  be  etched  in  the  backside  of  the  silicon  wafer  prior  to  the 
release  of  the  magnetometers.  These  lines  would  then  serve  as  crack  propagation  guides 
and  allow  for  each  individual  die  to  be  released  without  harming  the  magnetometers. 

The  first  prototype  PZT  xylophone  bar  magnetometer  devices  were  fabricated 
with  a  multilayer  (SiCVSijNVSiC^)  dielectric  stack.  The  goal  was  to  produce  a  released 
beam  with  zero  static  deflection.  These  devices  were  fabricated  but  still  experienced  a 
slight  degree  of  static  deflection  due  to  residual  stress  variability  in  the  multilayer  stack. 
However,  the  beams  could  still  be  tested  with  laser  Doppler  vibrometry  to  measure  the 
electrically  excited  (IV)  resonance  modes.  Figure  48  illustrates  a  few  of  the  resonance 
characteristics  for  one  of  the  devices  operating  at  atmospheric  pressure.  The  first  two 
resonance  modes  represent  a  spurious  non-symmetric  resonance  mode  that  was 
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associated  with  the  extremely  large  static  deflections  in  the  beam  ends.  The  fundamental 
resonance  mode  was  at  69.5  kHz  with  a  1st  harmonic  mode  at  162.1  kHz. 


Figure  48  Laser  Doppler  vibrometry  measurements  of  a  xylophone 
magnetometer  structure  of  the  4  lowest  resonance  peaks  (IV  excitation  at  atmospheric 
pressure).  The  top  images  correspond  to  asymmetric  resonance  modes  due  to  the 
deformed  ends;  while  the  bottom  2  images,  c  &d,  correspond  to  the  fundamental  and 
1st  harmonic  resonance  modes,  respectively. 


Figure  49  gives  the  frequency  scan  of  a  similar  device  (designated  pzt4)  tested  in 
the  APL  beam  deflection  rig  using  a  10  mV  excitation.  The  deflection  produced  by  this 
excitation,  applied  with  the  top  and  bottom  electrodes  connected  together,  was 
independent  of  applied  magnetic  field  and  in  close  agreement  with  the  resonance 
frequencies  obtained  with  the  laser  vibrometer.  The  piezoelectric  excitation  was  a  result 
of  the  different  resistances  of  the  top  and  bottom  electrodes  producing  a  net  voltage  drop 
across  the  PZT  film. 

A  high-resolution  scan  of  the  fundamental  resonance  peak  of  pzt  4  is  given  in 
figure  50.  The  Q-factor  is  very  low  at  113.  These  results  demonstrate  that  current  PZT 
devices,  because  of  their  large  piezoelectric  actuation,  are  unsuitable  for  detecting  static 
magnetic  fields. 
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Figure  49  Frequency  scan  (10  mV  excitation  at  l(f6  Torr)  of pzt4  device  in  APL 

beam  deflection  characterization  rig 


a. 


Figure  50  High-resolution  scan  of fundamental  resonance  peal  of  sample  pzt4. 


7. 


CONCLUSIONS 


* 


* 


The  objective  of  this  project  was  to  develop  a  sensitive  MEMS-based  magnetometer  for  use 
as  an  orientation  sensor  on  spinning  projectiles.  The  magnetometer  was  derived  from  the  JHU/APL 
resonating  xylophone  bar  with  sensitivity  linearly  proportional  to  the  drive  current,  mechanical  Q  at 
resonance,  and  the  magnetic  field.  Polysilicon  magnetometers  demonstrated  the  size  effects 
expected  for  such  devices.  Their  sensitivity  was  limited  by  the  relatively  high  sheet  resistivity  that 
restricted  the  drive  current.  A  CMOS-based  capacitive  sense  chip  was  developed,  fabricated  through 
MOSIS  and  tested.  This  circuit  can  be  used  in  future  designs  of  magnetometers  and  other  structures 
such  as  accelerometers  and  gyroscopes.  PZT-based  magnetometers  (processed  at  ARL,  Adelphi) 
were  shown,  at  the  current  state  of  development,  to  be  unsuitable  for  use  as  vibrating  bar 
magnetometers  because  of  low  mechanical  quality  factors,  asymmetric  conducting  layers,  and  high 
residual  stress  levels.  The  most  promising  technology  appears  to  be  based  on  CMOS-MEMS, 
although  more  development  is  required  to  produce  practical  devices.  Even  more  recent  advances  in 
silicon-on-sapphire  (SOS)  technologies  may  make  it  possible  to  replace  standard  capacitive 
transduction  schemes  with  much  more  sensitive  optical  ones.  Studies  have  shown  that  a  MEMS- 
based  magnetometer  would  have  important  applications  in  addition  to  orientation  sensing, 
particularly  when  combined  with  accelerometers  and  gyroscopes. 
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